Charge-ice dynamics in the negative thermal expansion material Cd(CN)2 
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We use variable-temperature (150-300 K) single-crystal X-ray diffraction to re-examine the in- 
terplay between structure and dynamics in the ambient phase of the isotropic negative thermal 
expansion (NTE) material Cd(CN)2. We find strong experimental evidence for the existence of 
low-energy vibrational modes that involve off-centering of Cd 2+ ions. These modes have the ef- 
fect of increasing network packing density — suggesting a mechanism for NTE that is different to 
the generally-accepted picture of correlated Cd(C/N)4 rotation modes. Strong local correlations 
in the displacement directions of neighbouring cadmium centres are evident in the existence of 
highly-structured diffuse scattering in the experimental X-ray diffraction patterns. Monte Carlo 
simulations suggest these patterns might be interpreted in terms of a basic set of 'ice-rules' that 
establish a mapping between the dynamics of Cd(CN)2 and proton ordering in cubic ice VII. 

PACS numbers: 65.40.De, 61.43.-j, 64.60.Cn, 77.84.-s 



Of the few materials known to exhibit the rare prop- 
erty of isotropic negative thermal expansion (NTE; vol- 
ume contraction on heating), cubic Cd(CN)2 shows the 
strongest such effect!^] Its volume coefficient of ther- 
mal expansion a = dlnV/dT = —61.2(12) MK" 1 is 
substantially more negative than the thermal expansion 
coefficients of typical engineering materials are positive 
(a ~ 30 MK" 1 )^ The origin of NTE in this material 
has never been identified conclusively, despite the re- 
lated material Zn(CN) 2 (which also shows NTE, but at 
a reduced level) having been the focus of a significant 
number of experimental and computational studiesP^^ 
Both compounds share the anticuprite structure in which 
tetrahedral Zn or Cd centres are co nnected via linear 
M-C-N-M linkages [Fig. [jja)].™^ Density functional 
theory (DFT) calculations suggest that the phonon rela- 
tions are similar in the two materials, with the more pro- 
nounced NTE effect of Cd(CN)2 a result of an enhanced 
pressure sensitivity of the low-energy mode frequencies.^ 
The atomic displacements associated with these crucial 
low-frequency modes are generally assumed to involve 
concerted rotations of M(C/N)4 coordination tetrahedra, 
which result in transverse displac ements o f the C and N 
atoms away from the M. . .M axesP^^llSJ Thermal pop- 
ulation of these modes would give rise to larger transverse 
displacements — hence reducing both the M. . .M separa- 
tion and the unit cell dimensions. 

Yet there are reasons to suspect that the dynamical 
behaviour in Cd(CN)2 and Zn(CN)2 might differ more 
strongly than is presently appreciated. It is known, 
for example, that Cd(CN)2 exhibits a structural phase 
transition at 150 K that is not observed in Zn(CN) 2 
at temperatures down to 20 K (albeit that the precise 
nature of this transition has not yet been reported) P 
This means that the Pn3m anticuprite structure used 
in DFT calculations-^ does not actually reflect the ther- 
modynamic ground state at K. And whereas tetrahe- 
dral geometry is the most common c oordin ation mo- 
tif for Zn 2+ ions in cyanide complexes j 13 * 14 ! the larger 
Cd 2 + cation (r ~ 0.95 A for Cd 2 + vs 0.60 A for Zn 2+ )±S] 
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FIG. 1: (a) The PnZm anticuprite structure of Cd(CN) 2 : 
Cd centres are shown as large red spheres, and cyanide ions 
as small black spheres. The cyanide orientations are disor- 
dered; consequently no distinction is made between C and 
N atoms. Like cuprite itself, two topologically independent 
networks interpenetrate. Here one is shown in ball-and-stick 
representation and the other in stick representation, (b) Cad- 
mium off-centering in Cd(CN)2. Displacement of the Cd 2+ 
ions parallel to one of the (100) crystal axes gives rise to a 
pseudo-octahedral coordination environment. 



nearly always adopts an octahedral coordination envi- 
ronment in otherwise equivalent compounds! 11 * 16 * ^ This 
difference in geometry preference is relevant because it 
has the potential to influence the relative importance 
of geometry-preserving 'rigid-unit modes' (RUMs) — 
precisely the type of ph onon m ode usually implicated in 
NTE for these materials 

In this study, we use variable-temperature (150-300 K) 
single-crystal X-ray diffraction to re-examine the inter- 
play between structure and dynamics in the ambient- 
temperature phase of Cd(CN) 2 . We find strong exper- 
imental evidence for low-energy vibrational modes that 
involve off-centering of Cd 2+ ions. These modes have the 
effect of increasing network packing density — so provid- 
ing an alternative mechanism for NTE. Strong local cor- 
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relations in the displacement directions of neighbouring 
cadmium centres are evident in the existence of highly- 
structured diffuse scattering in the experimental X-ray 
diffraction patterns. With the aid of Monte Carlo simu- 
lations we are able to interpret these patterns in terms of 
a basic set of 'ice-rules' that establish a mapping between 
the dynamics of Cd(CN) 2 and proton ordering in cubic 
ice VII. 

Our study makes use of a single crystal of Cd(CN) 2 
that was prepared by slow evaporation of a concentrated 
aqueous solution (see SI for details). X-ray diffraction 
patterns were collected using an Oxford Diffraction (Ag- 
ilent) SuperNovae diffractometer covering a full sphere of 
reciprocal space, with the same experimental procedure 
repeated at 30 K intervals between 150 and 300 K; refine- 
ments were carried out with CRYSTALS^ (see SI for 
further experimental details). While highly-structured 
diffuse scattering was clearly evident in each of these six 
data sets, the observed Bragg reflection conditions were 
consistent with the Pn3m structure reported in previous 
studies.^ On cooling below 150 K the bulk of this diffuse 
scattering intensity transferred into a large number of 
additional superlattice Bragg reflections (although some 
diffuse scattering remained) . The positions of these addi- 
tional reflections suggested an initial doubling, then qua- 
drupling of the unit cell in one or more directions. The 
complexity of this diffraction pattern and the likelihood 
of complex twinning associated with the phase transition 
means that we have restricted the scope of this particular 
study to the higher-symmetry ambient phase. 

The cubic unit cell parameters extracted from our 
diffraction data reflect the previously-reported NTE ef- 
fect [Fig. [2] P There are relatively few structural parame- 
ters to be refined crystallographically; a summary of the 
various values we obtain across the entire temperature 
series is given in Table SI. Here we focus on one crucial 
parameter: namely, the isotropic displacement parame- 
ter of the Cd atom, U (Cd). This parameter describes the 
mean-squared displacement of the Cd atom away from 
its crystallographic site, and is well constrained by X- 
ray diffraction intensities by virtue of the large Z con- 
trast between Cd and C/N atomsP^ Counterintuitively, 
we find that its magnitude actually decreases with in- 
creasing temperature [Fig. [2] . Working within the high- 
temperature approximation, we have 

U(Cd) = (u 2 ) = £g- t (1) 

where m is the mass of Cd and hcu c ff an effective (average) 
phonon energy. Consequently, a negative value of dU/dT 
implies a large and positive value of dw c ff/dT — i.e., the 
existence of strongly anharmonic low-frequency phonon 
modes, with the corresponding mode eigenvectors involv- 
ing Cd displacements. Such behaviour is commonplace in 
NTE materials, with the corresponding atomic displace- 
ment patterns generally in terpre ted as the microscopic 
driving force for NTE itself PBESSo the fact that Cd dis- 
placements have a strong and anomalous temperature de- 
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FIG. 2: Temperature-dependent structural evolution in 
Cd(CN)2 as determined using single-crystal X-ray diffraction. 
With increasing temperature, both the cubic lattice parame- 
ter a (open black circles) and the Cd mean-squared displace- 
ment [/(Cd) (filled red squares) decrease. The corresponding 
linear fits to data are included as guides-to-the-eye; estimated 
standard errors are smaller than the data points. 



pendence in Cd(CN) 2 suggests that these displacements 
are crucial to the NTE mechanism. We note that the ro- 
tational RUMs discussed in the introductory paragraph, 
and cited elsew here as a possible driving force for NTE 
in this materialpE^l do not give rise to Cd displacements. 

In order to investigate more fully the form of these dis- 
placement modes, we calculated difference Fourier maps 
for the Cd site for each of the single crystal X-ray diffrac- 
tion data sets at hand; the corresponding functions are 
illustrated in Fig. [3} At 300 K the modulat ions in the dif- 
ference function are small, indicating that there are no 
especially important regions of the unit cell where elec- 
tron density has neither been sufficiently accounted for 
nor overaccounted. For data sets collected at successively 
lower temperatures, a set of 'lobes' of unaccounted-for 
electron density emerges, with these lobes positioned at 
approximately I A away from the Cd site in each of the 
six symmetry-equivalent (100) directions. The picture 
that emerges resembles closely the behaviour observed in 
archetypal dynamic order-disorder transitions, such as 
the displacement of Pb atoms along (100) axes in cubic 
PbTi0 3 P 

Here in Cd(CN) 2 , displacements along the (100) direc- 
tions have an intriguing effect on the coordination geom- 
etry of the Cd 2+ ion. The displacement is towards an 
edge of the Cd(C/N)4 coordination tetrahedron, bring- 
ing the Cd centre closer to two of its four cyanide ligands 
(we refer to these as 'proximal' CN _ ions), and further 
away from two others [Fig. [ljb)]. In addition, a new as- 
sociation is now made with a further two cyanide ions 
belonging to the adjacent Cd(CN) 2 network (we refer to 
these as 'distal' CN~ ions). Consequently, the Cd off- 
centering we observe reflects a drive towards increased 
Cd coordination number: i.e. from four to six. At the 
same time, the cyanide ions undergo an effective increase 
in coordination number from two towards three as they 
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FIG. 3: Thermal evolution of residual scattering density 
around the Cd site in the difference Fourier map for Cd(CN)2. 
Each cube corresponds to an octant of the unit cell < 
(x,y,z) < |, with the Q, j, |) Cd position at its centre and 
the (100) crystal axes parallel to the cube edges. Isosurfaces 
are drawn at a level of p = +0.2 e~ A -3 . 



are now associated with an additional Cd 2+ in a sideways 
geometry [Fig. [ljb)]. The strong X-ray scattering con- 
trast between Cd and C/N atoms means that we are not 
especially sensitive to C and N displacements; however we 
note that there is clear evidence for transverse displace- 
ments of the distal cyanide ions towards the approaching 
Cd 2+ centre in the difference Fourier maps (see SI), sug- 
gesting possible anharmonic coupling of Cd translations 
along (100) directions with transverse CN displacements. 

Because thermal excitation of these modes increases 
the packing density (the unperturbed Pnim geometry 
being maximally extended) , each mode will be associated 
with a negative Griineisen parameter 7 = d In uj/d In V < 
— hence a driving force for NTEP Indeed a similar dis- 
placement pattern is observed in reverse as a mechanism 
for framework expansion in rutile-Ti02 under negative 
pressure: the (6,3)-net expands at — lOGPa to give a 
lower-density distorted tetrahedral TiC>2 net with the 
anticuprite topology.^ Distortive modes have also been 
suggested as a possible partial mechanism for NTE in 
the related framework NMe4[CuZn(CN)4], although in 
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FIG. 4: (Left) Experimental X-ray diffuse scattering patterns 
collected at 270 K for the (hkQ) and (hkk) cuts of reciprocal 
space and (right) the corresponding predicted diffuse scatter- 
ing from the Monte Carlo model as described in the text. 



that system the displacements involved are much smaller 
(<0.05A)P 

A displacement of ca 1 A is sufficiently large that one 
might anticipate strong correlation in displacement di- 
rections for neighbouring Cd atoms. Indeed, the exis- 
tence of such local correlations is supported experimen- 
tally by the observation of highly structured X-ray dif- 
fuse scattering [Fig. [4] and S1-S3]. The intensity of this 
diffuse scattering is temperature-dependent (see SI), in- 
dicating that its origin is dynamic rather than either 
static or compositional. 2 ® The scattering is also trans- 
verse polarised in that the intensity of the observed dif- 
fuse streaking is strongest along directions of reciprocal 
space perpendicular to the streaking itself; this is a sec- 
ond indication that its origin is displacive! 28 * 30 ! In gen- 
eral terms the scattering takes the form of diffuse planes, 
running across the (hkl)* ± e(h + l,h + l, 21)* regions of 
reciprocal space. It also obeys the extinction conditions 
h + k + l = 2n, allowing us to infer that the atoms respon- 
sible are separated by ^ (111) vectors in real-spacepS i.e., 
the Cd. . .Cd separation. For completeness, we note that 
while there are likely to be correlations between cyanide 
ion displacements, the large contrast in Cd vs C/N scat- 
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tering by X-rays will mean that diffuse scattering due 
to Cd displacements alone dominates the experimental 
data. 

So what is the simplest possible ordering model for 
Cd displacements that could explain the general form of 
this diffuse scattering? Our starting point in answering 
this question was the single assumption that the close 
approach between a given Cd 2+ centre and its proximal 
cyanides [Fig. [ija)] means that any cyanide ion is likely 
to be proximal only to a single Cd 2+ centre. In this 
way, the corresponding Cd-CN bond might relax to ac- 
commodate the initially-decreased Cd. . .C/N separation. 
It transpires that there exists a one-to-one mapping be- 
tween this rule and the problem of proton ordering in 
cubic ice: just as each O atom in the latter is associated 
more closely with two of its four H neighbours, so too is 
each Cd centre closer to two of its four cyanide neigh- 
bours in Cd(CN)2- And, just as each proton in cubic ice 
is covalently bonded to a single O atom, each cyanide 
ion in Cd(CN)2 is proximal only to a single Cd 2+ centre. 
Consequently, if this rule were to hold then for any given 
pair of connected Cd atoms their displacements must be 
'ice- rules' allowed [Fig. [5]. 

Making use of a simple Monte Carlo (MC) routine, 
we generated Cd(CN)2 atomistic configurations with 
Cd centres displaced according to this rule and subse- 
quently calculated the corresponding single-crystal dif- 
fuse scattering patterns with the approach described in 
Refs. [29f30l (see SI for details). We find that— on its 
own — this ice-rules-only algorithm predicts modulated 
diffuse scattering that is not yet convincingly close to the 
observed patterns [Fig. S5]. Consequently, our next step 
was to increase the complexity of the MC algorithm by 
considering explicitly the correlations between the direc- 
tions of nearest-neighbour Cd displacements. There are 
two types of Cd. . .Cd pairs separated by \ (111) vectors: 
'intraframework' pairs connected via cyanide ions, and 
'interframework' pairs which are not directly connected 
at all. By calculating diffuse scattering patterns for a 
range of different MC coupling constants J ln ter, -/intra we 
found the strongest resemblance to the experimental data 
for a model where intra-framework pairs displace in a 
like sense, and inter-framework pairs displace in opposite 
senses. Two representative diffuse scattering patterns are 
shown in Fig.[4j and full details of the various models are 
given as SI. 

The dominant Cd displacement mode is consequently 
an antifcrroclcctric (zone boundary) mode, correspond- 
ing to a local symmetry lowering from the space group 
Pn3m to IAi/amd: the displacements within either one 
of the two interpenetrated networks tend to occur in sim- 
ilar directions, but the direction is reversed from one net- 
work to the other. The distorted IA\/amd (\/2 x \[2 x 2) 
cell may well provide a useful starting point for structure 
solution of the low-temperature phase. Intriguingly, the 
same spa ce-g roup relation is found for the ice VII/ VIII 
transitionPB While the role of static disorder in c ubic i ce 
VII remains a controversial and difficult problem] 32 ! 33 ! it 
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FIG. 5: Correlated Cd 2+ displacements in Cd(CN) 2 . A re- 
quirement that each CN~ be proximal to a single Cd 2+ cen- 
tre implies ice-rules-like correlations in Cd 2+ displacements. 
Here, a displacement of one reference Cd 2+ centre (high- 
lighted) in a downwards direction reduces the number of pos- 
sible displacement directions from six to three for the Cd 2+ 
centre to which it is immediately connected within its cova- 
lent network. Of these three ice-rules-allowed directions, one 
is parallel to the original displacement vector (ferroelectric 
displacements), and two are perpendicular. Within MC simu- 
lations, the coupling parameter Jmtra can be used to favour ei- 
ther type of correlated displacement. Similarly, inter-network 
correlations can be controlled through the coupling parame- 
ter Jinter- The calculated diffuse scattering patterns shown in 
Fig. [4] were generated using a pair of Jintcr, Jintra values that 
favour ferroelectric displacements within the same network, 
but antiferroelectric displacements between networks. 



is pertinent to note that the phase is thought to exhibit 
isotropic NTE for some regions of its stability fieldP^ As 
for ice itself, one anticipates that Cd(CN)2 must exhibit 
a large entropic signature associated with Cd displace- 
ment ordering at the 150 K phase transition. We have 
not attempted to measure the temperature dependence 
of the specific heat as part of this study, noting that the 
added complication of CN orientational disorder may ren- 
der specific heat anomalies difficult to interpret.^ 

Returning to the displacement pattern shown in 
Fig. [ljb), we note that the antiferroelectric modes we 
have identified will mean that pairs of unconnected Cd 2+ 
centres displace in opposite senses, giving rise to localised 
polar dimers (effectively a Cd2 quadrupole). That is, 
the proximal cyanides of one Cd 2+ ion become the distal 
cyanides of the corresponding inter-framework Cd 2+ ion 
and vice versa. Very similar behaviour has been observed 
in frustrated niobate pyrochlores such as Y2NbTi07 and 
Y2Nb207, where correlated Nb off-centering is thought 
to be driven by the propensity of Nb 4+ to form bonded 
metal-metal dimers ES1311 Whereas the Nb-Nb contacts 



5 



are fixed and hence the systems represent a static 'charge 
ice', in cubic Cd(CN)2 the charge-ice displacements are 
clearly dynamic. 

In summary, we have shown that Cd off-centering 
is a dominant dynamical process in the NTE material 
Cd(CN)2- Displacements of the Cd 2+ ions away from 
their high-symmetry (crystallographic) site have the ef- 
fect of increasing their coordination number and also re- 
sult in a contraction of the cubic unit cell. We find 
that experimental X-ray diffraction patterns for Cd(CN)2 
contain highly-structured diffuse scattering. This diffuse 
scattering, which has a dynamical origin, can be inter- 
preted in terms of correlated displacements of Cd 2+ cen- 
tres according to a set of local antiferroelectric 'ice-rules'. 
The existence of charge-ice dynamics in Cd(CN) 2 lends 



further weight to the emerging correspondence between 
geometric frustration and NTE! 22 * 38 l 39 lthe same configu- 
rational degeneracy responsible for entropic anomalies in 
cubic ice is also implicated in driving NTE in Cd(CN)2- 
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